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ABSTRACT 

We constrain blastwave parameters and the circumburst media of a subsample of ten Bep- 
poSAX Gamma-Ray Bursts. For this sample we derive the values of the injected electron energy 
distribution index, p, and the density structure index of the circumburst medium, k, from simul- 
taneous spectral fits to their X-ray, optical and nIR afterglow data. The spectral fits have been 
done in count space and include the effects of metallicity, and are compared with the previously 
reported optical and X-ray temporal behavior. Using the blastwave model and some assump- 
tions which include on-axis viewing and standard jet structure, constant blastwave energy and 
no evolution of the microphysical parameters, we find a mean value of p for the sample as a 
whole of 2.04±°;° 2 . A statistical analysis of the distribution demonstrates that the p values in 
this sample are inconsistent with a single universal value for p at the 3tx level or greater, which 
has significant implications for particle acceleration models. This approach provides us with a 
measured distribution of circumburst density structures rather than considering only the cases 
of k = (homogeneous) and k = 2 (wind- like). We find five GRBs for which k can be well 
constrained, and in four of these cases the circumburst medium is clearly wind-like. The fifth 
source has a value of < k < 1, consistent with a homogeneous circumburst medium. 

Subject headings: gamma rays: bursts 



Introduction 



Since the discovery of GRB afterglows in 1997 (jCosta et alJll997t IVan Paradiis et al.lll997l ; iFrail et al 



19971 ). their relatively longer duration and broader wavelength coverage compared with the prompt emis- 
sion has made them the most accessible and arguably the most profitable area for observational GRB 
studies. Their non-thermal spectra and lightcurves can generally be well described by t he fireball model 
(jCavallo fc Reesl Il978t iRees k Meszarosl Il992l; iMeszaros k Reesl Il997t IWiiers et all Il997t) which describes 
a decelerating, highly relativistic outflow, the so-called blastwave, interacting with a surro unding external 
medium. Adaptations hav e been made to accommodate a non-unifo rm external medium ( Meszaros et alj 



1998 



Chevalier fcllll999h and the fact that the outflow is collimated |Rhoadslll997l Il999[ ISari et al.lll999h 



The evidence for collimation comes from energy considerations, and directly from observations of an achro- 
matic break in afterglow lightcurves, the so-called jet-break. The jet-break is an indication of the time at 
which sideways spreading of the jet begins to become important, combined with the fact that the effect of 
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relativistic beaming starts to be noticeable within th e jet. For the l ong-duration GRBs, which most prob- 



ably originate in the cor e-collapse of a massive star (|Wooslevlll993t iMacFadven et al.ll2001t IWooslev et al 



20021; iHiorth et al.ll2003h . a n(r) oc r 2 density profile from a stellar wind would not be unexpected for the 
external medium. Its density can be probed through the temporal decays of the afterglow and, contrary 
to expectations given their stellar progenitors, long GRBs studied thus far do not all require a wind-like 
n(r) oc r~ 2 density profile in the local external medium (e.g. iPanaitescu fc Kumarll2002f ). 



The synchrotron nature of the blastwave radiation produces spectra and lightcurves comprising of a 
set of power laws with characteristic slopes and frequencies. Accurate measure ments of these observed 
quantities allow the underl ying parameters of the blastwave to be determined (e.g. Wiiers fc Galama 1999t 
Panaitescu fc Kumarll200lf ) . Optical through X-ray spectra provide the opportunity to measure the index of 
the input power law energy distribution of electrons, p, potentially the index k of the density profile of the 
circumburst medium, n(r) cx r~ k , and in some cases also the cooling break frequency v c — the frequency of 
electrons whose radiative cooling time is equal to the dynamical timescale of the blastwave. 

If the microphysics of all GRBs is assumed to be the same, the range of values of the input electron energy 
distribution indices should be narrow. However, the directly measurable data that lead to the parameter p, 
such as the break frequencies and power-law slopes of the spectra and lightcurves, are dependent on various 
other factors, like the circumburst density, the fraction of energy in electrons (e e ) and magnetic fields (es), 
and simply the total amount of energy, making it more difficult to obtain a consistent value for p. Further, 
determination of these parameters so far works most successfully for the later-time afterglow considered here: 
before approximately 0.1 days the observed lightcurves and spectra are likely some mixture of the prompt 



emissi on (a ttributed to internally colliding shocks) and afterglow (the external shock), see e.g.. lNousek et al 



( 2006 ) and O'Brien et al. ( 20061 ). These authors discuss the 'canonical' X-ray lightcurve for Swift bursts 



which has three phases, the final phase beginning between 0.01 and 0.12 days after trigger (jNousek et al 



2006) and showing the type of steady decay seen in the pre-Swift era. The advantage of using late-time data 
only is that we can reliably restrict ourselves to the slow cooling regime in modeling the blastwave, where 
the injection frequency of the electrons is well below thei r cooling frequency. The measurements we perform 
in this study are at late enough epochs (see Table 1 in Istarling et al. 2007 . hereafter Paper I), that they 
should not be affected by any prompt emission components. 

A potential problem is that occas ionally values of p below 2 have been found (e.g. for GRB 010222, 



Masetti et al.1 l2001[ IStanek et al.l l200lh . This requires a cut-off at the high-energy end of the distribution 



of th e electrons, and adaptations have been made for such cases (|Bhattacharvall2001t IPanaitescu fc Kumar 
20011 ). The details of these adaptations, however, are still under debate, since the evolution of the high- 
energy cutoff is not well constrained. Since the lowest values for p we find are ~ 2, we do not take any 
adaptations for this effect into account. 

For several bursts studied here the underlying param eter set has been measur e d independently. De- 
terminations for s e ts of GRBs have been made by e.g. IWiiers fc Galamal (|19991 ); IPanaitescu fc Kumar 
|200ll ); lYost et al.l |2003h : iGendre fc Bo^ (|2005h : iGranot et al.l J2OO6I). Unfortun ately, only a small frac- 



tion of GRBs have measurements in all wav ebands, most notably GRB 970508 (e.g. I Wiiers fc Galama 



199' 



Berger et al 



9; 



Panaitescu fc Kumar||2 002: Yos t et alJl2003L Van der Horst & Wijers in preparation) and GRB 030329 (e.g. 



2003bt IWillingale et al.ll2004 iFrail et al.ll2005l ). In the Swift era such studies are generally lim- 



ited to part of the parameter set, since they often use only one waveband, and therefore lack the possibility 
of finding the characteristic break frequencies in the broadband spectrum. 



Here, we fit the broadband spectral energy distributions (SEDs, from nIR through X-ray) of a subsamplc 
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of the BeppoSAX sample of GRB afterglows. We constrain a subset of the blastwave parameters, namely the 
index of the power law energy distribution of electrons, p, the density profile parameter of the circumburst 
medium, fc, and for some bursts the position of the cooling break, v c . Because of the paucity of radio data 
for most bursts in our sample, we have not included these, ensuring a more homogeneous approach. Radio 
data are not needed to determine p and fc, since this is possible from only the nIR through X-ray SEDs. 

We make use of simultaneous fits in count space to obtain the most accurate possible measurements. 
In Paper I we provide details of the observations, data reduction and fitting method, summarized here in 
Sections 2 and 3 together with the description of the model used. In Section 4 we present the results of our 
p- and fc-value and cooling break analysis, both for the sample and for individual sources. We compare these 
results to those of previous studies of this kind in Section 5, and perform statistical modeling of the p- value 
dataset. Here we discuss our findings in the context of the fireball model and long GRB progenitors. We 
conclude by summarizing our results in Section 6. 



2. Observations and Spectral Fitting Method 



X-ray observations were made with the narrow field instruments on-board BeppoSAX (Paper I, Table 
2), and here we have combined data from the MECS2 and MECS3 instruments (except in the case of 
GRB 970228, where we use the MECS3 instrument only). 

Optical and nIR photometric points were taken from the literature (Paper I, Table 3) and from our own 
observations of GRB 990510 (Curran et al. in preparation). Tempo ral decay s l opes w ere again taken from the 
literature ( see Paper I, Tab l e 1): the optical temporal slop e s from 



Zeh et al 



(2006) and the X-ray temporal 



slopes from lGendre fc Boerl (|2005l ) : lin 't Zand et alj (|l998h ; lNicastro et al] ( 19991 ) . We have transformed the 
time of the optical and nIR photometry to the log of the midpoint of the combined X-ray observations. We 
avoid using data taken before 0.1 days after the GRB, hence we assume no complex and flaring behavior 
occurs in the lightcurves and we restrict ourselves to the slow cooling regime. All data are transformed to 
count space for fitting purposes, in order that no model need be assumed a priori for the X-ray spectrum, 
and fitted within the ISIS spectral fitting package ijHouck fc Denicolall2000l ). Models consist of either a single 
or a broken power law, to allow for a possible cooling break in between the optical and X-rays. In the broken 
power law model we fix the difference in spectral slope to 0.5, as expected in the case of such a cooling 
break. Both Galactic and intrinsic absorption are components in the models, allowing for either Milky Way, 
Large Magellanic Cloud (LMC) or Small Magellanic Cloud (SMC) extinction laws for the GRB host galaxy, 
at either Solar (Z ), LMC (#=1/3 Z ) or SMC (Z=l/8 Z ) metallicity. All details of the observations, 
reduction and analysis are given in Paper I, together with the results of the power law plus extinction fits. 



3. Theoretical Modeling 

We assume that the ambient medium density as a function of radius can be described as a power law with 
index fc, i.e. n(r) oc r~ k , so that a homogeneous medium is given by fc = and a stellar wind environment 
by fc = 2 — the two most likely scenarios. For the purpose of looking at optical and X-ray data at ~ 1 day, 
we need to derive the time dependency of the peak flux F Vt7nax and the cooling frequency v c as a function 
of fc, and the peak frequency v m (which has no dependence on fc), assuming the afterglows are in the slow 
cooling regime. After the jet break all of these parameters have no dependence on fc; in this region we know 
we are dealing with a jet geometry hence we label this case 'Jet' or j (Table [lj. 
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k n 1 
2(4-fe) u 2 


-1 
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4-3fe 1 1 
2(4-fc) 2 2 





v m 


3 3 3 
2 2 2 


-2 



Table 1: The temporal power law indices of the peak flux F u . maXl the cooling frequency v c and the peak 
frequency v m , as a function of the circumburst density profile index k for pre- (columns 2-4) and post- 
jet-break (column 5). 
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If one assumes that the flux is a power law in frequency and time with /3 (or F) the spectral slope and 
a the temporal slope, using the conventions F v oc v ~@t ~ a oc v ~( r ~ 1 H _Q , with power law photon index 
r = f3 + 1 , one can derive these slopes as a function of k and the power-law index p of the electron energy 
distribution. These values for a, (3 and T are given in Table [5] for two different situations: the observing 
frequency in between v m and v Cl and the observing frequency above both frequencies. Also shown in this 
table are the closure relations between a and 0. 

One can invert all these relations to obtain p from a and /3, and even determine k from these observables: 

= 4(3/? ~2a) = 4[3(r - 1) - 2a] 

3/3-2a-l 3r-2a-4 ' 1 j 
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Table 2: The temporal and spectral slopes of the flux, a and /3 (or T, where r = (3 + 1) respectively, the 
closure relations between a and /3 (or T), and p as a function of a, /3 and T. 
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From Table [2] it is clear that when the observing frequency is higher than v m and v c , both a and (3 only 
depend on p, and do not depend on k. In the situation where the observing frequency is situated in between 
v m and v c , the spectral slope only depends on p, but the temporal slope depends on both p and k. So the 
structure of the ambient medium can only be determined in the latter case {y m < v < v e ), although having 
more accurate information on p from the situation with v m < v c < v is useful to get a better handle on k. 

4. Results 

The results of fits to the SEDs for all GRBs in the sample are given in both Table 4 and Fig. 2 of Paper 
I. For derivation of the blastwave parameter p we adopt the best-fitting models as listed in Paper I and in 
particular cases additional models were included. SMC-like absorption was the preferred extinction model in 
all cases except for GRB 000926 where LMC-like absorption is preferred. We calculate the values for p and 
k for two cases: the cooling frequency in between the optical and X-rays, and the cooling frequency above 
both. We have checked whether the cooling frequency could lie below the optical band using the relations 
of the fireball model, but this turns out not to be the case for these GRBs. The cooling frequency, v c , is 
obtained directly from the SED fits for 3 bursts: GRBs 990123, 990510 and 010222, with v c = (1.3±^;|) x 10 17 , 
(4.3±0.5)xl0 15 , and (4.1^ 5 1 2 ) x 10 15 Hz at 1.245, 1.067 and 1.511 days since burst respectively. Applying 
the fireball model we find that another two sources, GRBs 980329 and 980703, may require a break within 
their SEDs, at - 2.6 x 10 17 and ~ 8 x 10 17 Hz at 1.148 and 1.333 day s since burst (the logarithmic midpoint 
of the X-ray spectrum) respectively. 

The resulting values for p can be found in Table [3] and for k in Table [4] All errors are given at the 
90% confidence level for one interesting parameter, unless otherwise stated. a± and a,i refer to the pre- and 
post-break optical lightcurve slopes given in Paper I; we allowed for the possibility that these breaks are not 
jet breaks by considering that oti is both pre- (columns 6 and 7) and post- (column 8) jet break. This has 
also been done for the X-ray temporal slopes. 

4.1. Results: Individual Sources 

4.1.1. GRB 970228 

We find p = 2.44t°;og and k = 1.73t°;gg, with ^ x < v„ at 0.52 days. A cooling break between 
optical and X-ray bands is not required at the time of our SED (0.52 days) according to the F-test: the 
F-test probability, the probability that the result is obtained by chance, is 2.1 x 10~ 2 , which is quite high; 
so adding one extra free parameter is not a significant improvement. Using the best fitting model of a 
single power law plus SMC-like extinction, in the regime i>x < v c , we find that the data can be fit by both a 
homogeneous and a wind-like circumburst medium; the value of k is best constrained by the optical temporal 
slope. 

4.1.2. GRB 970508 

We find p = 2.56^0 46 and k — 0.49^0 67) with < v c at 1.68 days. Using the single power law with 
LMC extinction (rather than SMC to obtain the more conservative errors on the spectral slope) and optical 
to X-ray offset free, in the regime i>x < we find the data are best fit by a homogeneous medium; the value 
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k is best constrained by the optical temporal slope. Broadband modeling by Van der Horst & Wijers (in 
preparation) with k as a free parameter gives very tight constraints on k: a value of 0.02 ± 0.03 is derived, 
i.e. a homogeneous medium. 

Previous works put the cooling br eak at optical frequenc ies, v c = 1.6 x 10 14 Hz, at 12.1 days since 
burst, between the B and V bands (e.g. Wiiers fc Galamalll999fh We, however, find that the cooling break is 



likely to lie above the X-rays at 1.68 days. We note in this context the uncertain extrapolation of the optical 
data used in the SED, owing to an irregular shaped lightcurve at early times, which we have attempted 
to account for in allowing the optical to X-ray offset to go free. This is a particularly difficult case given 
that the X-ray data cover the time period immediately following an optical flare when the optical lightcurve 
appears to have flattened before breaking to a typical and well defined power law. 



4.1.3. GRB 971214 



We find p = 2.20 ± 0.06 and k = 2.17; 



-0.25 



with vx < v c at 1.36 days. Using the best-fitting model 



of a power law plus SMC extinction, in the regime vx < v c , the data are best fit by a wind medium: from 
the optical temporal slope we find k = 2.17^0 35, whilst from the X-ray temporal slope k = 2.331q'34- A 
broken power law does not provide a significant improvement in the fit compared to a single power law, i.e. 
the F-test probability is high (4. 1 x 10~ 2 ). A spectral brea k is claimed for this burst in the IR (~1 micron) 
at 0.58 days (y c ~ 3 x 10 14 Hz, iRamaprakash et al.lll998| ). This is not the cooling break, but the peak of 
the SED moving to lower frequencies, so there is no conflict with our results for v c . 



4.I.4. GRB 980329 

We find p = 2.50t£e2 and k = -4.&9±¥?£ s , with v c < vx at 1.15 days. The spectral fit obtained with 
a single power law plus SMC extinction is inconsistent with the optical temporal slopes. A spectral break in 
the power law does not provide a significant improvement of the fit according to the F-test (probability of 
7.2 x 10~ 2 ), but this spectral break model provides agreement between the spectral slopes and the optical 
and X-ray temporal slopes. In this regime, v c < vx, k cannot be well constrained though the centroid of 
the fit to the optical data is that of a homogene ous medium. We note that when omitting the I band point 
from the SED, which may be overestimated (see Yost et al. 20021 and Paper I), our results do not change. 



4.1.5. GRB 980519 



We find p 



2.96±g;g| and k 



0.23; 



-1.22 

-3.05' 



with vx < v c at 0.93 days. Using the power law plus SMC 



extinction model for the SED, in the regime vx < v ci the optical data are best fit by a homogeneous medium, 
and the X-ray temporal slopes can be fitted by both a homogeneous and a wi nd medium; k is therefore best 
constrained by the optical temporal slope. In contrast, Chevalier fc Li ( 19991 ) found that the radio emission 
of the afterglow of GRB 980519 measured between 7. 2 hours and 63 da ys since the burst is consistent with 



an external wind instead of a homogeneous medium. Frail et al. (2000) note, however, that the interstellar 



scintillation present in the radio data does not allow one to draw firm conclusions on this. 

The optical temporal break at 0.48±0.03 days ( Zeh et alll2006[ ) cannot be explained by passage through 
the optical bands of v Cl since the derived value for p from the temporal slope is too high in that case 
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(p = 3.69 ± 0.06) compared to the p-value from the joint spectral slope. It also cannot be explained by a 
jet break, since p is too low in that case (p = 2.27 ± 0.05). It appears that the fireball model is a good 
explanation for the first temporal slope and the spectrum used here, but the post-break optical slope has 
either been incorrectly measured or we do not yet have the correct model for this afterglow. We note that 
this afterglow showed a very steep temporal decay compared to other GRBs (jHalpern et al.lll999h ■ This is 
somewhat reminisce nt of the very early-time decays of many Swift bursts which occur at <500 s after trigger 



( Nousek et al. 20061 ) and likely have a significant prompt emission component. 



4.1.6. GRB 980703 

We find p — 2.74lg'4g and k = 1.63^5 6 3 4 6 , with v c < vx at 1.33 days. The spectral fit obtained with 
a single power law plus SMC extinction is inconsistent with the X-ray temporal slopes. A spectral break in 
the power law does not provide a significant improvement of the fit according to the F-test (probability of 
7.2 x 10 -2 ), but this spectral break model provides agreement between the spectral slopes and the optical 
and X-ray temporal slopes. Hence we use the broken power law plus SMC extinction model in the regime 
Vc < Vx- The nature of the optical temporal break at 1.35 ± 0.94 days cannot be determined because of 
large uncertainties in the optical temporal slopes, which are also the reason why k cannot be constrained. 



Two publications have postulated a position for th e cooling bre a k in past studies. IVreeswijk et al 



(|1999| ) propose v a < vx < v c at 1.2 days after the burst. iBloom et al.l (|1998f) propose v Q < v c < vx at 5.3 
days, and our SED study at 1.33 days, when compared with the optical temporal slope, is inconclusive since 
both vx < v c and vx > v c can be accommodated. It may be that the cooling break has moved into our 
observed bands during accumulation of the X-ray spectrum (possibly indicated by the inability of the fireball 
model to fit the data when a single power law is assumed for the spectrum). If we require consistency with 
these previous results, the cooling break must be moving to lower energies and lies approximately at X-ray 
frequencies in our data. This would mean that the circumburst medium is homogeneous, since v c is expected 
to move as t -1 / 2 in this case, while v c will increase in time as i 1 / 2 in the wind case. 

The host galaxy of GRB 980703 appears to have a high and possibly variable optical extinction along 
the line of sight to the GRB (see Paper I). The different (and formally inconsistent) values of Ay may be 
due to different methods for measuring the extinction, probing of different regions of the host galaxy, or 
may indicate that the environment in which the burst occurred is changing with time. We have used the 
Vreeswiik et al. ( 19991 ) optical data and scaled it from 1.2 days to 1.33 days after trigger. Combining the 



op tical and X-ray data wh en fitting provides us with a different estimate for the extinction than was obtained 



by IVreeswiik et al.l (|1999l ) for the optical alone. Any change in optical extinction will have an effect upon 



the measured spectral slope and hence the derived value of p. 



4.1.7. GRB 990123 

We find p = 1.99±^ and k = 2.00±^;|f , with v c < vx at 1.25 days. With this best fitting SED model 
of a broken power law plus SMC extinction, in the regime v c < i>x, the optical temporal decay is best fit by 
a wind medium. The optical temporal break at 2.06 ± 0.83 days is marginally consistent with a jet-break: 
the p- value derived from the post-break temporal slope is consistent with the one derived from the pre-break 
optical temporal slope, but inconsistent with the spectral slope. The uncertainties in the X-ray temporal 
slope are too large to determine the phase of blastwave evolution (i.e. before or after jet-break) from the 
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X-ray data alone. 



For GRB 990123 the temporal slope difference between optical and X-ray of 0.25 also agrees with the 
spectral analysis where a broken power law model is the best fit, indicating a cooling break between the 
optical and X-ray bands at 1.25 days post-burst . The value we der ive for p is consistent with that derived 
from the X-ray spectrum alone of p — 2. ±0.1 (Stratta et al. 2004 ) and lower than a previous estimate via 
broadband modeling of p — 2.28 ± 0.05 ( Panaitescu fc Kumar 2002f) (we note that the latter uncertainty is 



1 a and not the 90% error used in the rest of our paper). 



GRB 990510 



We find p = 2.06: 



0.02 



and k 



0.80^q'qq, with v c < vx at 1.07 days. The best fitting model to the 



SED is clearly a broken power law with negligible extinction, in the regime v c < vx- There is considerable 
improvement in the x 2 when allowing for this break in the power law, noted by previous author s, which we 
find is located at 0.016 - 0.020 keV at ~1.07 days since burst (consistent with the value given bv lPian et al 



20011 ). The slope change is as expected for a cooling break in the slow cooling regime when leaving both 
power law slopes free. 

The optical data are best fit by the values for p and k mentioned above. In this case, however, the X-ray 
temporal slope is not consistent with the spectral slope nor the optical temporal slope at the 90% (~ 1.6cr) 
level, but does agree within 3er. The optical temporal break at 1.31 ±0.07 days is marginally consistent with 
a jet-break. The der ived value of p i s consi stent with the value derived from the BeppoSAX X-ray spectrum 
alone of p ~ 2.1, by Kuulkers et al. ( 2000l). and r ules out the value of p = 2.6 ± 0.2 also derived from the 
BeppoSAX X-ray spectrum bv lStratta et al.l (|2004 ): that same X-ray data is used here, but is combined with 
nIR and optical data to obtain our limits on p. 



4.1.9. GRB 000926 

We find p = 2.54+°;Jg and k = 2.16^%, with is x < v c at 2.23 days. With the power law plus LMC 
extinction model, in the regime vx< v c , the optical temporal decay is best fit by a wind- like medium. Large 
uncertainties in the X-ray temporal slope prevent determination of the circumburst medium structure and 
blastwave evolution phase from X-ray data alone; the optical temporal break at 2.10 ±0.15 days is consistent 
with a jet-break. 



4.1.10. GRB010222 

We do not find a consistent solution for this afterglow taking 90% uncertainties, but we do find one 
taking 3a uncertainties: p = 2.04+°;^ and k = 2.28^0 29, with v c < vx at 1.51 days. Adopting the single 
power law model with LMC extinction, in the regime vx < v c , the optical slopes are not consistent with 
the spectral slope nor the X-ray temporal slopes, both at the 90% and 3<r levels. Using a broken power 
law plus SMC extinction, in the regime v c < vx, for which the F-test indicates a marginal improvement 
(probability of 1 x 10~ 4 ), the X-ray temporal slope is not consistent with the spectral slope nor the optical 
temporal slope at the 90% level, but they are consistent at 3<r. In the first case we obtain p = 2.721"q'q5 and 
k = 0.30lQg4; in the latter case p = 2.041° and k = — 2.25l4 3 54 , which is derived from the pre-break 



2.28^0 29, derived from the post-bre ak optical tempo ral slope 



optical temporal slope, or k 

temporal break happens quite early, 0.64 ± 0.09 days after the burst (|Zeh et al 



Since the 
200fih . and the post-break 



optical slopes are inconsistent with a jet-interpretation, the early temporal slope is probably influenced by 
late-time energy injection and a medium with k — 2.28^Q2g the correct interpretation. 

Panaitescu fc Kumar ( 2002h find the cooling break to lie at optical wavelengths or longer, in agreement 



with our results. However, they derive a low value for p of 1.35, and find significant reddenin g of the optical 



spectr u m which they say explains the second steepening observed in the optical after 6 days bv lFruchter et al 



(|200l[ ). iBhattacharva fc Resmil (|2004f ) obtain good fits by adopting a two-slope electron energy distribution 
with p- values of 1.3 and 2.1, below and above a so-called injection break, respectively. This injection break 
is located in the X-ray regime at ~ 1 days after the bur st. In their model the optical temporal break is 
a jet-break and the circumburst medium is homogeneous. iBiornsson et al.l (|2002l ) argue that the unusually 
slow decay of this afterglow and positive detection of linear polarization can be explained by a jet model 
with continuous en ergy injection. Such slow decays are seen in the 'canonical' Swift X-ray lightcurves (e.g. 



Nousek et al.l 120061 ) . termed the plateau phase, and typically begin a few hundred seconds after the GRB 
trigger. The electron energy distribution then has p = 2.49 ± 0.05, which is inconsistent with our result for 
p from the spectral fits. 



5. Discussion 



The parameters that can be derived from broadband modeling of afterglow lightcurves describe the 
micro- and macrophysics of the relativistic jet and its surrounding medium. To obtain the full set of 
parameters the spectral energy distribution has to be covered from X-ray to radio wavelengths. Two of these 
parameters can be deduced from just the spectral and temporal slopes in the optical and X-rays, i.e. the 
electron energy distribution index p and the circumburst medium profile parameter k. These two parameters 
have been determined in this paper for a selection of 10 GRBs, for which the final results are shown in Table 
[5] and Fig. Q] For completeness the values for v c are also listed in Table El half of the GRBs in this sample 
have a value v c < vx ~ 8 x 10 17 Hz. 



Table 3. Values for p. We calculate the results for the cases k = and k = 2. Bold type highlights consistent results (at the 90% level). 
In cases where the best-fitting spectral model to the SED (from Paper I) is inconsistent with the model fits, we show the results for this 

original best-fitting model in italics. 
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Table 4: Values for fc. Bold type highlights consistent results (at the 90% level). In cases where the best- 
fitting spectral model to the SED (from Paper I) is inconsistent with the model fits, we show the results for 
this original best-fitting model in italics. 
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Fig. 1. — Derived values of p (top panel) and k (lower panel) for each individual afterglow (see Table [5]): 
the horizontal axes represent the GRBs in date order left to right and errors are 90 % confidence. In the top 
panel the dotted line indicates the most likely value of p = 2.04; in the lower panel the dotted lines indicate 
k = (homogeneous medium) and k = 2 (stellar wind). 
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Table 5: Final results for p, k and v c for all ten bursts in the sample. 



"Optical data extrapolation is uncertain 

6 Broken power law better than single power law in SED. 

c This solution is consistent with all measurements except the post-break optical temporal slope. 
d Brokcn power law better than single power law in SED. Large uncertainties in the optical temporal slopes. 
e X-ray temporal slope only consistent at 3cr level. 
^X-ray temporal slopes have large uncertainties. 

9 Break in the optical lightcurvc at 0.6 days is not a jct-brcak. X-ray temporal slope only consistent at 3c level. 
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5.1. The Distribution of p 



Some theoretical studies of partic le acceleration by ultra-relativistic shocks indicate that there is a nearly 



5 

universal value of p of 2.2 — 2.3 (e.g. iKirk et al.ll2000t lAchterberg et al 



that there is a large range of possible values for p of 1.5 — 4 (e.g. lBarinsll2004r ). From the results presented 



2001 ). while other studies indicate 



in this paper and from broadband modeling of individual bursts, quite a large range of values for p have 
been found, which could indicate that there is a large intrinsic scatter in the value of p. Here we test the 
null- hypothesis namely that the observed distribution of p can be obtained from a parent distribution with 
a single central value of p, by performing a statistical log-likelihood analysis on the obtained values of p. 

We first determine the most likely value of p for our sample, by minimizing the log-likelihood of our 
10 measured values of p. In order to do this, we describe the measured values of p and their uncertainties 
as probability distributions, in which we take asymmetric measurement uncertainties into account. These 
probability distributions are given as two halves of normal distributions with different widths that are set 
by the measured uncertainties. The two halves of the probability distributions are normalized such that 
they are connected as a continuous function and have a total integral equal to 1 . This results in distribution 
functions which are given by 



V2 



V^ci + <J 2 ) 



3 -(p-Pmca B ) 2 /2< 
= -(p-Pmca B ) 2 /2<T2 



(P < Pn 
(P > Pn 



(2) 



with iTi and 02 the lower and upper 1 a uncertainties in the measured value of p (indicated as p mea s) respec- 
tively. This probability function describes how likely a value of p is given the measurement (p m eas,f 1,02)- 
To convert the 90% confidence limits in Table [5] to la uncertainties, we divided those by a factor of 1.6. 
The log-likelihood for these probability distributions is given by 



JY 
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(3) 



in which in the last term is used for p < p me as,z and a^.i forp > p mea s,i- N is the number of measurements, 
i.e. in this case N = 10. 



The log-likelihood is a function of p and we determine the most probable p-value by minimizing this 
function for our 10 bursts. We find the most likely value of p = 2.04+QQg, with a log-likelihood of 613.6. 
We have derived the uncertainties in this most likely p-value by generating 10 5 synthetic datasets for the 
10 bursts in our sample. These datasets are possible measurements within the measurement uncertainties, 
generated by taking random numbers from the probability distributions that are defined by the measured 
p- values and their uncertainties. For each possible dataset we determine the most likely p- value, and from 
the resulting distribution of most likely p-values we obtain the 1 a uncertainties in the value of p = 2.04. 

To test our null- hypothesis that the observed distribution of p can be obtained from a parent distribution 
with a single central value of p, we generate 10 5 different synthetic datasets for p for the 10 bursts in our 
sample by taking random numbers from probability distributions that are described by the most likely value 
of p = 2.04 and the 1 a uncertainties in the 10 values of p. We then take 10 values of p coming from these 
synthetic datasets (one random number for each measurement), calculate the most likely value of p and the 
accompanying log-likelihood. The resulting distribution of log-likelihood values is plotted in Fig. [21 together 
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with the minimal log-likelihood of 613.6 coming from the measurements. We find that in 99.92 % of the cases 
the log-likelihood of the synthetic data is smaller than the one coming from the measured p-values. This 
leads to the conclusion that the hypothesis, that the distribution of p from our sample can be obtained from 
a parent distribution with a single central value of p, is rejected at the 3 a '-level. 

This result challe nges the theoretical studies of particle ac celeration which claim that there is one uni- 
versal valu e of p (e.g. iKirk et al.l 12000 : lAchterberg et al.l 120011) . and is consistent with similar findings by 
Shen et al.l (|2006l) from fits to t he prompt emission of a sample of Burst And Transient Source Experiment 
(BATSE) GRBs. iBarina (|2004T ) shows that this often quoted value of p = 2.23 is a special case, with partic- 
ular assumptions on hydrodynamic quantities, like the compression ratio of the upstream and downstream 
velocities, on the influence of the magnetic field on the dynamics, and on electron scattering angles. He 
claims that there is a large range of possible values for p of 1.5 — 4. We derive constraints on the width of 
the distribution of p-values below. 
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Fig. 2. — The log-likelihood distribution coming from the 10 5 synthetic datasets generated from one single 
value of p (solid line); the dashed vertical line indicates the log-likelihood for the real data. In 99.92% of 
the cases the log-likelihood of the synthetic data is smaller than the one coming from our measured sample. 
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From Table [5] and Fig. Q] it can be seen that the log-likelihood is dominated by the value for p of 
GRB 990123, and after that by GRB 980519. In fact, just the two p- values of GRB 990123 and GRB 980519 
are different by 12 a, and thus they alone refute the constant p hypothesis. We also examined the strength 
of the evidence against a universal p from samples excluding one of the 10 bursts in our sample. Only in 
cases where either GRB 990123 or GRB 980519 are excluded is the significance less than 3er, but even then 
the hypothesis can be rejected at the ~ 2 c-level. 

Since we now know that the measured values of p from our sample are not coming from a parent 
distribution with a single central value of p, we can put constraints on the width of the parent distribution of 
p- values. We do this by introducing an intrinsic scatter on the most likely value of p = 2.04: we generate 10 5 
synthetic datasets by first taking random numbers from normal distributions with a central value of p = 2.04 
and a certain scattering cr SC at; we then take random numbers from the probability distributions that are 
centered at these scattered p-values and have a width given by the measurement uncertainties. Performing 
the same log-likelihood analysis as described above for varying <r sca t , we find that the 3 a result holds for 
Cscat > 0.03; at the lcr level this is cr sca t > 0.24. We can also constrain the upper limit on the scatter in 
the parent distribution of p-values from our sample. This can be derived by determining for which cr sca t 
the log-likelihood of the measured p-values is smaller than the one from the synthetic dataset analysis in 
99.73 % of the 10 5 cases. This upper limit for tJ sca t is 1.40 at the 3 a level; at the 1 a level this is <x sca t < 0.45. 
In summary, from our sample we constrain the intrinsic width of the parent distribution of p-values to be 
0.24 < cr scat < 0.45 at the 1 a level and 0.03 < cr scat < 1.40 at the 3 a level 

Values of p — 2.2 have been widely quoted as a typical number and the deviations from this interprete d 
as due to the externa l environment or further energy injections from the source (e.g. iBerger et al.ll2003ah . 
Meszaros et al. I (|l998h showed that for a single value of p, variations in the jet energy per solid angle, i.e. 
structured jets, could lead to a range of lightcurve decays; w hich can also be pr oduced if a set of intrinsically 
similar, structured jets are viewed at various off-axis angles (jRossi et alJl2002l ). The study of such structure 
in the jet and viewing angle dependence is, however, beyond the scope of this work. 



Johannesson et al.l (|2006l ) claim that from broken power law fits the value of p is overestimated from the 
pre-break lightcurve slope while being underestimated from the post-break slope, if the temporal coverage is 
limited. The results presented in that paper were only for a homogeneous ambient medium; they claim that 
in a wind medium the p value from the pre-break lightcurve slope does not show a systematic deviation. 

More recently, breaks in lightcurves, both optical and X-ray, which would generally have been attributed 
to jet-breaks, have been found to be chromatic, in sharp contrast to the picture of a jet-break which should 
be achromatic. As such, several previously claimed jet-breaks may have to be revised; unfortunately, there 
is no relevant data to confirm or reject these claimed jet-breaks, since these have been found only in one 
waveband (mostly the optical, owing to the lack of dense monitoring in X-rays in the pre-Swift era). As 
such, the mention of jet-breaks in this paper has been taken at face value, but with the caveat mentioned 
here. 

If we compare the values for p calculated here with those calculated from the X-ray spectra alone 



(jStratta et al.l 120041 ) . we find they agree within the 90% errors except for the burs ts GRB 970508 and 



GRB 990510. With this method we reduce the average 90% error on p from ±0.58 (jStratta et all 120041 ) 



+0.12 
-0.20 



, and the values obtained arc also likely to be more robust given that consistency between the nIR, 



to 

optical, UV and X-ray is required. 



Shen et al. ( 20061 ) have measured the minimum p ossible width of the electr on energy index distribution 
for X-ray afterglows of BeppoSAX bursts taken from |Pe Pasauale et al.l (|2006l ). by fitting single absorbed 
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power laws to the spectra and deriving p. They found that the narrowest possible distribution is consistent 
with a delta function within the la errors, in contradiction to their findings from a Swift sample of X- 
ray afterglows and from the prompt emission of a BATSE GRB sample. They comment that the relatively 
larger errors on the X-ray slopes of the BeppoSAX sample compared with the Swift sample are likely allowing 
for an apparently narrower distributio n. They calculate that the distribution of p for the Sw i ft sam ple of 
X-ray afterglows (jO'Brien et al.l 120061 ) has a width of 0.34 ± 0.07. We note that Ishen et all |2006h have 
derived distributions of p from X-ray afterglows only, which makes it impossible to know whether the cooling 
frequency lies above or below the X-ray regime, and can only be determined if there are accompanying 
optical measurements. 

Panaitescu fc Kumar ( 2002h calculated jet parameters for a sample of 10 GRBs including several studied 



here, using broadband observations including radio data when possible. In their model, based also on the 
fireball model, they assume uniform jets (no structure across the jet) and the the energy parameters e e and 
€b are constant, and finally they assume the observer is located on the jet symmetry axis. Our spectrally- 
derived values are consistent with theirs at the 3a level or better for GRBs 970508, 980519, 990510 and 
000926. There is no agreement for GRB 990123. They found a spread in p values as do we, but with half 
the values lying below 2, and a mean value of p = 1. 92^20 (2ct). 



Chevalier fc Lil (|2000h carried out a similar study of broadband afterglow data, and their estimates for p 
are in agreement with ours for the GRBs 970228, 970508, 980519 and 990510, and disagree for GRB 990123 
(they do not quote errors per source but estimate errors to be ^0.1). They conclude that their sample shows 
a range in the values of p which is not consistent with a single value. 

We note that in all these studies the results are dependent upon the chosen model and simplifying 
assumptions. 



5.2. The Circumburst Medium 



The profile of the circumburst medium is a particularly important parameter in studying the progenitors 
of GRBs. In the case of long-soft bursts the progenitor is a massive star that is expected to have had a 
massive stellar wind in earlier phases of its evolution, and one would expect to see a signature of that wind 
in the afterglow lightcurves. Evidence for a stellar wind in the form of fast outflowing absorption lines within 
restframe UV spect r a has been seen in a small number of cases, the best example being GRB 021004 (e.g. 
Schaefer et al.ll2003l : IStarling et al.1 12005'). This does not mean, however, that a density profile with n cx r~ 2 
is always expected, since this assumes a constant mass-loss rate and a constant wind velocity. Changes in 
mass-loss rate and also interactions of the wind with the interstellar medium can alter the ambient medium 



profile (e.g. iRamirez-Ruiz et al.ll2005l : Ivan Marie et al.ll2006l ). 



In the first broadband modeling attempts the ambient medium was assumed to be the homogeneous 
interstellar medium, which was consistent with the derived particle densities. However, since the progenitors 
of at least a fraction of all GRBs are now known to be massive stars and the blastwave is situated at ~ 10 17 
cm at ~ 1 day after the burst, a massive stellar wind profile is expected. Nonetheless, the medium that the 
blastwave is probing could still be homogeneous, because of the emergence of a reverse sh ock in the wind when 
the wind meets the interstellar medium (see e.g. Wijersl 2001 : Ramirez-Ruiz et al. 2005 ). This shocked wind 
turns out to become homogeneous and, for typical physical parameters derived from afterglow modeling, the 
blastwa ve encounters the tran sition from the stellar wind to this homogeneous shocked medium at ~ 1 day 
(see e.g. lPe'er fc Wiiersl koOfi). The actual time of the transition, which would be detectable in the afterglow 



-21 - 



lightcurves, depends for instance on the mass-loss rate and the density of the interstellar medium, which are 
both not really well constrained for most GRBs. In our sample we do not see this kind of transition in the 
optical lightcurves in which there is a break. It has been claimed for some bursts, for instance GRB 030329 



(jPe'er fc Wiierd 120061 ), that this transition is observed, but the often complex structure of the lightcurves 



confuses the determination of such a transition. 

Another way to obtain a constant density from a massive stellar wind is in the region after the wind 
termination shock. The distance to the termination shock can be very large, but recent observations of two 
Wolf-Rayet bin aries has suggested that this distance could be several times smaller if the wind is asymmetric. 



aiyi 
i'e ( 



Eldridgei (|2007f ) shows that wind asymmetry probably exists in two systems, that can be caused for example 
by rotation, which is expected for GRB progenitors in the framework of the collapsar model in order to retain 
enough angular momentum. If the asymmetry exists for the entire stellar lifetime, then a closer termination 
shock and asymmetric supernova may be expected, increasing the chances of observing an afterglow traversing 
a constant density medium. 

In our sample of 10 GRBs there are four sources that are consistent with an r~ 2 wind medium, with 
relatively small uncertainties, namely GRBs 971214, 990123, 000926, and 010222. There are three GRBs 
which arc not consistent at a 90% confidence level with a wind medium, GRBs 970508, 980519, and 990510, 
although for GRB 970508 caution is warranted with the interpretation of the lightcurves; and for the other 
four bursts the uncertaintie s are t oo large to constrain the ambient medium profile. We contrast our findings 



with iPanaitescu fc Kumar! (J2002J), who, in broadband fits to the data of 10 bursts, found that a wind-like 
medium was preferred over a homogeneous medium in only one case: GRB 970508. For this particular burst 
our analysis provides a value of p which is an equally good fit to wind or uniform medium predictions for 
v < v c from the spectra and lightcurves, but the closure relations are obeyed (at the 2a level for both 
cases) if v > v c . They find circumburst densities of order 0.1-100 cm~ 3 for most sources, which they argue 
demonstrates the surrounding medium does not have, in general, the r~ 2 profile expected for the unperturbed 
wind of a massive GRB progenitor. 

The associati on of long-soft GRBs with Ib/c supernovae was first suggested for GRB 9 80425 / SN 1998bw 
by Galama et al. |l998h . and confirmed by the discovery of GRB 030329 / SN 2003dh (e.g. lHiorth et al.ll2 003). 



Therefore, it is useful to compare the circumburst medium characteristics derived from GRB afterglows and 
from radio observations of supernovae, which also trace the density profile of the surroundings of these massive 
stars. Around radio supernovae r~ 2 density profiles have been found, but also in some cases significantly 



flatte r behavior of ~ r 15 in SN 1993J and SN 1979C (for a review on radio supernovae see IWeiler et al 



2004 and references therein for individual supernovae). In the latter case a transition from r to r 
was even observed. This flatter density profile can be attributed to changes in the mass-loss rate of the 
massive star in some phases of its evolution. The three sources in our sample of GRBs with a density profile 
flatter than r~ 2 are possible examples of the relativistic blastwave ploughing its way through a region of the 
circumburst medium which is affected by changing mass-loss rates. So although in Table[S]we have described 
them as GRBs with a homogeneous ambient medium, this is not necessarily the case. Especially for GRBs 
970508 and 980519 this is a possibility, but the uncertainties on k are too large to distinguish a homogeneous 
from a flattened wind medium. GRB 990510 has smaller uncertainties and seems to be more consistent with 
a homogeneous medium, especially since the upper limit on k is ~ 1.0, which is much flatter than what is 
observed in radio supernovae. 

The Swift satellite now provides us with substantially greater coverage of a large number of X-ray 
afterglows (^100 per year) and often with high quality data from which to measure the spectral and temporal 
slopes. However, few of these also have substantial optical follow-up. The combination of X-ray and optical 
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data helps determine the position of the cooling break and obtain accurate spectral slopes which provide the 
value of p. For the derivation of k in this study, we have found the optical temporal data most constraining. 
For this reason, and for the confirmation of achromatic jet breaks it is essential that such late-time optical 
data be obtained for as many well sampled Swift X-ray afterglows as possible. 

6. Conclusions 

We have measured the injected electron energy distribution index, p, in the framework of the blastwave 
model with some assumptions which include on-axis viewing and standard jet structure, constant blastwave 
energy and no evolution of the microphysical parameters. We have also measured the density profile of the 
circumburst medium, n(r) cx r~ fc , from simultaneous spectral fits to the X-ray, optical and nIR afterglow 
data of 10 BeppoSAX GRB afterglows. 

A statistical analysis demonstrates that the distribution of p values in this sample is inconsistent with 
a single value for p at the 3 a certainty, which is at variance with many theoretical studies of relativistic 
particle acceleration. We constrain the width of the parent distribution of p values and find it to be of the 
order of a few tenths, with p — 2.04^'^ as the most likely p-value in our sample. 

We measure the distribution of the local density parameter k, generally only assumed to be or 2, 
and we find that the majority of GRBs for which we can constrain k well are consistent with a wind-like 
circumburst medium. One source (GRB 990510) is clearly, i.e. > 3tr certainty, inconsistent with this picture 
and fits instead a homogeneous medium. These results are consistent with the expectations of at least a 
subset of GRBs originating from massive stars, which have a lot of mass-loss in the form of a surrounding 
stellar wind. We have discussed the possibility of values of k < 2 within the stellar wind framework. 

The method presented here provides a way to study the distribution of blastwave parameters for a 
sample of GRBs, and allows estimates to be derived when insufficient data are available for a full time- 
dependent solution. In the current Swift era the method is equally well applicable, although one would have 
to ensure that the data are in the afterglow domain, i.e. not contaminated by prompt emission or late-time 
energy-injection. A decent sampling of the optical SED and lightcurve, more difficult with the average fainter 
Swift afterglow sample, is crucial to constrain the temporal slopes and cooling break frequency, which in 
turn are necessary to determine p and k. 
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